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OETECTION  PERFORMANCE  OF  NORMALIZER  FOR  MULTIPLE  SIGNALS  SUBJECT 


TO  PARTIALLY  CORRELATED  FADING  WITH  CHI-SQUARED  STATISTICS 


INTRODUCTION 


In  a  recent  study  [1],  the  detection  performance  capability  of  a 
multiple-pulse  system  subject  to  correlated  fading  was  quantitatively 
delineated.  It  was  assumed  there  that  the  noise  level  was  known,  so  that  a 
threshold  could  be  set  for  an  arbitrarily  specified  false  alarm 
probability.  Then  the  detection  probability  was  evaluated  as  a  function  of 
the  threshold  level,  the  received  signal-to-noise  ratio,  the  number  K  of 
signal  pulses,  and  the  fading  statistics. 

Here  we  will  extend  these  earlier  results  to  cover  the  case  where, 
additionally,  the  noise  level  is  unknown  and  must  be  estimated  on  the  basis 
of  a  finite  number  L  of  noise-only  samples.  The  same  approximation 
technique  that  was  presented  in  [1]  is  used  to  determine  the  detection 
probability  of  this  normalizer  system.  The  reader  is  referred  to  [1]  for 
additional  background,  motivation,  interpretations ,  and  related  references. 
For  the  sake  of  brevity,  we  will  employ  the  same  notation  and  presume  that 
the  reader  has  complete  familiarity  with  the  earlier  material  and 


deve 1 opment . 
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PROBLEM  DEFINITION 

We  will  couch  the  problem  in  a  particular  setting,  one  with  obvious 
appeal  and  application;  however,  it  should  be  obvious  how  to  extend  this 
setting  to  a  more  general  one,  particularly  in  light  of  the  arbitrary  fading 
covariance  coefficients  that  are  allowed  in  the  analysis. 

Suppose  a  sequence  of  K  tone  bursts  at  a  common  center  frequency  are 
transmitted,  as  depicted  in  figure  1.  Each  rectangular  slot  symbolizes 

-fV«Huencjj 


Figure  1.  Time-Frequency  Occupancy  Diagram 

a  tone  of  duration  T^  seconds  and  approximate  frequency  bandwidth  1/T^ 

Hz.  These  bursts  may  be  abutting  in  time  or  may  be  arbitrarily  separated  in 
time  by  several  multiples  of  T  .  At  the  receiver,  K  narrowband  filters  of 
bandwidth  1 /T ^  Hz  are  sampled  at  the  times  of  peak  signal  output  (if 
present)  and  their  squared  envelopes  are  summed.  Depending  on  the  time 
separation  between  pulses,  the  signal  strength  may  fade  considerably;  the 
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exact  amount  and  frequency  of  the  fading  depends  on  the  distribution  of  the 
fading  and  the  covariance  of  the  fading  amplitude  of  adjacent  (as  well  as 
separated)  pulses. 

It  is  presumed  that,  during  a  single  tone  of  duration  seconds,  the 
fading  is  essentially  constant,  resulting  in  a  constant  amplitude  scaling 
and  phase  shift  applied  to  the  pulse.  The  time  separations  between  pulses 
in  figure  1  are  arbitrary,  thereby  allowing  for  an  arbitrary  degree  of 
correlation  between  the  fading  factors  applied  to  each  pulse. 

To  establish  a  reference  against  which  this  sum  of  K  matched  filter 
outputs  can  be  compared,  for  purposes  of  deciding  on  the  presence  or  absence 
of  signal,  a  group  of  L  nonoverlapping  noise-only  slots,  located  arbitrarily 
in  the  time-frequency  plane,  are  also  energy-detected  and  summed.  For  very 
large  I,  this  noise  reference  is  very  stable,  and  performance  approaches 
that  predicted  by  [1].  However,  for  moderate  values  of  L  and  for  small 
false  alarm  probabi 1 ities  of  interest,  it  is  important  to  know  how  much 
degradation  in  performance  is  incurred  by  being  forced  to  use  this  noisy 
ref erence . 

An  obvious  implementation  of  the  processing  implied  by  figure  1  is  to 
employ  fast  Fourier  transforms.  The  L  reference  bins  can  then  be  an 
arbitrary  collection  of  time  and/or  frequency  bins.  However,  L  cannot  be  so 
large  that  nonstationary  and/or  nonwhite  noises  cause  their  own  kind  of 
errors  in  noise  power  estimation.  The  tradeoff  between  these  conflicting 
requirements  will  be  assessed  quant i ta t i ve 1 y  in  this  investigation. 
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DEFINITIONS  OF  PARAMETERS 

Very  heavy  reliance  will  be  made  here  on  the  basis  that  was  set  up  in 
[1].  Thus  we  have  the  following  fundamental  parameters  of  the  detection 
procedure  (the  immediate  references  In  tables  1  and  2  are  to  [1]): 

K,  number  of  potential-signal  pulses  added,  (figure  1  and  A-ll ) ; 

m,  signal  fading  parameter  (power-scaling  is  chi-squared  with  2m 
degrees  of  freedom),  (13); 

jp^,  normalized  covariance  coefficients  of  signal  power-scalings 

hi-  <I5>- 

average  received  signal  energy  per  pulse _ ,  (9); 

single-sided  received  noise  spectral  density  level 

number  of  noise-only  pulses  added. 

Table  1.  Fundamental  Parameters 


In  addition,  there  are  two  very  useful  auxiliary  parameters  that  found 


frequent  use  in  [1  ] : 


K 


e’ 


kj 


=  equivalent  number  of  independent  signal  pulses,  (10); 


N,  m  Kg  =  a  summary  parameter  describing  the  distribution  of 
the  sum  of  power  scalings, (A-24)  and  (A-29). 


Table  2.  Auxiliary  Parameters 


None  of  the  parameters,  m,  Kg,  N,  need  be  integer.  Also,  N  can  be  larger 
or  smaller  than  K,  the  number  of  signal  pulses. 

PROBAB 1 1  IT IES  FOR  KNOWN  NOISE  LEVEL 


The  probability  density  function  of  the 
signal  envelope-squared  samples  is  given  by 


sum  y  [1 ;  ( A-l 1 )  ]  of  the  K 
[1;  (B-4)  ] 


Py(u) 


exp(-u/a)  u 


K-l 


aK'N  bNr«) 


ifi(n:  “(1  -  0) 


for  u  >  0  , 


(1) 


where  [1 ;  ( A-32) , ( B-3) , ( B-7 ) ] 


‘-2%-  b  -  ?«n  ('  -  R>  •  '-ri 

o 


(2) 


The  exceedance  distribution  function  0  ( u )  of  output  sum  y  is  given  by 

Y 

several  alternative  forms  in  (1;  ( B -9 ) , ( B - 1 1 ) ,  (  B-l  3) ] .  For  a  fixed 
threshold  (known  noise  level),  the  detection  probability  is 
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PQ  =  (y.A,R,N,K)  = 


=  1 


1 


(1  +  R) 


I  ^  W  D 

n=0 


-  E 


K-l  4-n 


«]  ■ 


JL= 


2a. 


and  the  false  alarm  probability  is  [1;  (B-10)] 

PF  =  EK-1^  ’ 

where  we  define  the  exceedance  distribution  function 

En(x)  =  exp(-x)  en(x)  , 


(3) 


(4) 


(5) 


en(x)  =  ^  xj/j!  (6) 

3=0 

is  the  partial  exponential  [2;  6.5.11].  The  results  in  (3)  and  (4)  should 
be  used  for  L.  -  «>,  that  is,  for  known  noise  level. 


NORMAI I ZER  RATIO 

From  this  point  on,  L  is  presumed  finite.  Suppose  a  noise  level 
estimate,  y  ,  is  obtained,  based  upon  L  independent  measurements  of 
noise-only  bins.  It  is  assumed  that  the  average  noise  level  in  these  L  bins 
is  the  same  as  in  the  K  potential -signal  bins,  but  that  this  noise  level  is 
unknown.  If  we  let 
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Y  =  Y(K,E1/No) 


(7) 


denote  the  sum  of  K  signal  bin  outputs  with  average  signal -to-noise  ratio 


E../N  ,  then 
1  o 


Y0  =  y(L,0) 


(8) 


is  the  corresponding  sum  of  L  noise-only  bins.  Now  define  ratio 


V  Yq  Y(L,0) 


(9) 


for  sets  of  K  and  L  pulses,  respectively.  The  noise  contributions  to  the 
total  of  K  +■  L  outputs  are  presumed  independent  of  each  other;  however,  the 
signal  fading  factors  amongst  the  K  signal  outputs  are  correlated  to  an 
arbitrary  degree.  We  are  interested  in  the  distribution  of  this  normalizer 
ratio,  v. 


When  signal  is  absent,  the  ratio  u  in  (9)  is  independent  of  the  absolute 
level  of  the  received  noise;  therefore,  we  can  expect  the  normalizer  to 
achieve  the  important  capability  of  constant  false  alarm  rate.  That  means 
that  a  specified  false  alarm  probability  can  be  achieved  without  knowledge 
of  the  average  noise  power  level. 


The  quantities  y  and  yq  are  the  sums  of  K  and  L  squared-envelope 


samples,  respectively,  and  are  not  the  averages  of  these  sampled 
quantities.  In  terms  of  the  sample-average  quantities,  we  could  define  a 
slightly  different  normalizer  ratio 
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V 


L 

K 


V  . 


(10) 


It  then  readily  follows  that  the  cumulative  distribution  function  of  random 
variable  v  is 


Pp(u)  =  Prob(t>  < 


u)  =  Prob^u  <  i  oj  =  ?v(^  u)  , 


(11) 


in  terms  of  the  cumulative  distribution  function  of  ratio  v  in  (9).  Thus,  a 
simple  scale  factor  change  allows  for  consideration  of  the  alternative  ratio 

V  . 


When  we  plot  the  detection  probability  versus  the  false  alarm 
probability,  that  is,  eliminate  the  threshold,  the  same  performance 
characteristics  result  for  random  variable  v  as  for  Accordingly,  we  will 
not  use  or  refer  to  ~  or  P_(u)  any  further,  but  concentrate  solely  on 
normalizer  ratio  v,  given  by  (9). 


NORMAL  I  Z  E.  R  DISTRIBUTIONS 


The  characteristic  function  of  noise-only  random  variable  yQ  can  be 
found  directly  from  [1;  ( A - 1 3 ) ]  by  setting  A  to  zero  and  replacing  K  by  L: 


f  (?)  =  d  -it  a) 
'o 


(12) 
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The  corresponding  probability  density  function  of  y0  is 


p  (u)  ,  a!'1  »p(-u/a)  for  u  ,  0 
r0  T(L)  aL 


The  exceedance  distribution  function  is 


Q  (u)  =  Prob(Y  >  u)  =  £  (u/a)  for  u  >  0  , 

'o  L 


in  terms  of  the  functions  defined  in  (5)  and  (6). 


The  cumulative  distribution  function  of  ratio  v  in  (9)  is  given  by 
(since  yq  >  0) 

P^(u)  =  Prob(u  <  u)  =  ProbA*  <  =  Prober  <  uyq^  = 


=  j  dy  Pv(y)  dx  p  (x)  =  l  dy  p  (y)  Q  (y/u)  = 
J  «  J  'rx  *  T  T* 


for  threshold  u  >  0,  where  we  used  (1)  and  (14).  We  now  expand  E 
according  to  (5)  and  (6)  and  integrate  term-by-term, to  obtain  [3;  7.621  4] 


P  (u)  = 


(*f  (rh)K  3.  jr^  f ("•  K  :  ('  -  0  Hh:)  <  '6: 


But  f rom  (2) , 


miw,ttr.TC-w*'  ,~V-  ~ fi' 


where  the  parameters  involved  are  described  in  tables  1  and  2.  Making  these 
substitutions  in  (16),  there,  fol  lows  for  the  cumulative  distribution  function 
of  random  variable  u. 


P  (u) 

V 


K  ;  K; 


R  u 
1  +  R  1  -t-  u 


An  alternative  more  useful  form  is  obtained  when  we  use  [2;  15.3.3]: 


(19) 


for  u  >  0.  This  result  is  very  attractive  since  the  negative  integer 
argument ,  K  in  the  hypergeometric  function  causes  termination  of  the 
series  at  X  terms.  Thus,  (19)  is  a  closed  form  (albeit  tedious)  for  the 
cumulative  distribution  function  of  v,  involving  a  finite  number  of 
elementary  functions. 


2 

It  should  be  noticed  that  the  absolute  noise  level  a  does  not 

n 

appear  in  (18)  or  (19).  (The  cumulative  distribution  function  for 
alternative  normalizer  ratio  u  given  by  (10)  can  now  easily  be  found  by  use 
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COMPARISON  WITH  EARLIER  RESULTS 


The  result  (19)  for  the  cumulative  distribution  function  of  normalizer 
ratio  v,  operating  in  a  partially  correlated  fading  environment,  is  an 
approximation,  having  been  based  upon  a  characteristic  function  fitting 
procedure  explained  in  [1;  ( A-24) -(A-28) ] .  Nevertheless,  (19)  is  identical 
with  the  exact  fading  result  for  a  related  normalizer  problem;  namely, 
agreement  with  [4;  (25)]  is  achieved  under  the  following  identifications: 


TR  4783  Here 

a  U 

M  K 

N  L 

u  +•  1  N 


R 


Interpretation 

threshold 

number  of  signal  pulses 
number  of  noise-only  pulses 
m  Ke,  table  2 


Table  3.  Identification  of  Variables 


The  identity  of  v  +  1  with  N  is  made  by  comparing  [4;  ( 24 A )  ]  with 
[1;  '  A -29 ) ] .  The  final  identity  of  »  with  R  utilizes  [4;  (248)]  and 
[1;  (9)]: 


Rt  Et/N  E,  K/N 

T  x  T _ o  1 _ o 

w  V  +  1  N  N 


R  . 


where  the  arrow  indicates  transferrance  from  [4]  to  [1], 


(20) 
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The  approach  in  [4]  proceeded  as  follows:  the  detection  probability 
for  nonfading  signals  in  all  the  bins  depended  only  on  the  total  received 
signal-to-noise  ratio  R.J..  When  RT  was  assigned  the  fading  probability 
density  function  [4;  (24A)],  the  average  detection  probability  in  [4;  (25)] 
resulted.  For  the  special  case  of  fading  parameter  v  =  M  -  1  there, 
numerous  graphical  results  were  given  in  [4;  figures  1-36]. 

The  current  results  here  are  more  general,  in  that  they  allow  for 
partially  correlated  fading  (through  parameter  Kg)  and  a  more  general 
power-fading  model  (with  2m  degrees  of  freedom).  This  means  that  N  =  m  Kg 
here  is  not  restricted  to  be  equal  to  the  number  of  signal  pulses,  K,  but  is 
arbitrary.  Thus  the  current  numerical  results  will  significantly  augment 
and  extend  those  in  [4].  If  N  *  K  here,  then  R  =  TTJ/No  =  signal-to-noise 
ratio  per  pulse,  and  (19)  reduces  to  [4;  ( 1 5B) ] ,  for  which  many  numerical 
results  were  given  in  [4;  figures  1-36], 


SPECIAL  CASES 


For  m  =  1,  which  corresponds  to  Rayleigh  amplitude  fading,  and  for 

p  =  A  ,  which  corresponds  to  uncorrelated  fading,  then  K  =  K,  N  =  K,  and 
r  J  K  J  6 

we  get  from  ( 19) , 


in  agreement  with  [4;  ( 1 5 B ) ] . 


(21 


On  the  other  hand,  if  R  -  0,  then  (18)  and  (19)  both  reduce  to 
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L-l 


P(0)(u) 

V 


-  W 


(Kk 


X'.  (1  +  uj* 


(22) 


which  is  equal  to  1  -  Pp,  where  Pp  is  the  false  alarm  probability. 


2  . 


Since  noise  level  on  is  not  involved  in  (22),  threshold  u  can  be 


selected  to  realize  a  given  Pp,  once  K  and  L  have  been  specified.  This  is 
a  quantitative  verification  of  the  expected  constant  false  alarm  rate 
property  of  the  norma lizer. 


Finally,  in  the  special  case  of  one  signal  pulse,  K  =  1,  and  Rayleigh 
amplitude  fading,  m  =  1,  then  Kg  =  1 ,  N  =  1 ,  R  =  E^/N^,  and  (19)  yields 

L-l 


.<»>  ■  rrvrr  JL  (ufhj  - 1  -  (rl  l  * »)' 


(23) 


Jl= 0 


That  is, 


-  ■  (  ♦  rh)  '  • 


1 


(24) 


which  agrees  with  [5;  (6)]  when  we  make  the  identifications  (from  there  to 
here)  of  N  -»  L,  T/N  -»  u,  y  -»  R . 


RECURSION  FOR  CUMULATIVE  DISTRIBUTION  FUNCTION 


let  the  hypergeometric  function  appearing  in  (19)  be  represented  as 
fol lows : 


(K) 


H|(x)s  F(-Jl,  K  -  N;  K;  x)  . 


(25) 
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^H^(x)  =  [K  +  21  -  2  +(N  -  K  +  1  -i)x]  H  ^(x)  -  (K  +  JL  -  2)(1  -  x)H^_2(x) 

for  £>  1  ,  (27) 

where  we  define  H  (x)  =  0.  In  terms  of  (25),  the  cumulative  distribution 
function  of  v  in  (19)  becomes 


This  form,  in  conjunction  with  recursion  (27),  was  used  for  all  the 
numerical  results  here,  for  L  finite.  The  parameters  appearing  in  (28)  have 
all  been  explained  in  tables  1  and  2.  The  explicit  dependence  on  the 
fundamental  parameters  is  indicated  below: 

Ke  =  Ke(K,  fpkj})  , 

N  =  N(m,  K,  fpkj})  . 

R  =  R(E7/N0,  m,  K,  ?pkj})  .  (29) 

In  addition,  the  cumulative  distribution  function  in  (28)  is  a  function  of  L 
and  threshold  u. 


14 


DETECTION  ANO  FALSE  ALARM  PROBABILITIES 

The  detection  probability  is  given  by 

PQ  =  Prob( v  >  u J R  >  0)  =  1  -  Pu(u)  ,  (30) 

where  P„(u)  is  available  in  (28).  The  false  alarm  probability  is 

Pp  =  Prob( v  >  u)r  =  0)  =  1  -  P(v0)(u)  ,  (31) 

where  P^°^(u)  is  available  in  (22).  By  allowing  threshold  u  to  vary 
over  a  wide  range,  PQ  and  Pp  values  can  be  obtained  and  plotted  against 
each  other,  resulting  in  the  standard  receiver  operating  characteristics; 
the  threshold  is  thereby  eliminated  from  the  plotted  outputs.  Programs  for 
plotting  PQ  vs  Pp,  both  for  L  finite  as  well  as  infinite,  are  listed  in 


appendix  A. 
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graphical  results 


Due  to  the  multitude  of  parameters  appearing  in  this  investigation  (see 
tables  1  and  2),  it  is  impossible  to  give  a  comprehensive  compilation  of 
encompassing  numerical  results.  Considering  just  the  covariance 
coefficients  {pkjl -|  for  the  moment,  complete  specification  requires 
assignment  of  K(K  -  l)/2  values  to  these  quantities;  to  circumvent  this 
difficulty,  we  consider  numerically,  here,  only  the  very  special  case  of 
exponential  correlation,  for  which 


for  1  <  k,  j  <  K 


(3 


and  look  at  a  couple  of  particular  values  for  p.  Our  approach  here,  of 
necessity,  is  to  give  some  representati ve  sample  receiver  operating 
characteristics  and  a  general  computer  program  in  BASIC,  whereby  additional 
results  can  easily  by  obtained  once  the  user  has  specified  all  the 
particular  values  of  interest  in  his  application.  This  program  allows  for 
arbitrary  covariance  coefficients,  [p^  ,  and  is  not  limited  to  the 
spec  i  f ic  example  ( 32) . 


The  particular  cases  we  will  investigate  are  as  follows: 

K  «  1.  2.  4, 

L  =  16,  32,  «0, 


TR  8133 


V . 


All  possible  combinations  of  these  four  fundamental  variables  lead  to  30 

•k 

plots,  which  appear  below  in  figures  2-31.  (There  are  only  6  plots  for 
K  =  1,  not  12,  because  the  value  of  p  is  irrelevant  for  K  =  1).  The  curves 
are  indexed  by  the  per-pulse  signal-to-noise  ratio,  r^N  ,  in  dB.  The 
false  alarm  and  detection  probability  pairs  range  from  (poor  quality)  pair 
(.5,.oi)up  to  (high  quality)  pairs  near  (l  E-l  0,  .  999). 

The  number  of  signal  pulses,  K,  is  limited  to  the  low  values  1,  2,  4, 
because  these  seem  to  be  the  cases  of  most  immediate  practical  use.  The 
number  of  noise-only  samples,  L,  is  not  evaluated  for  L  =  64  because  of  the 
proximity  of  the  results  to  those  for  L  =««;  conversely,  results  are  not 
presented  for  L  =  8,  because  a  severe  degradation  in  performance  occurs,  that 
probably  cannot  be  tolerated.  The  fading  parameter  value  m  =  1  corresponds 
to  Rayleigh  amplitude  fading  (exponential  power  fading),  while  m  -  .5 
corresponds  to  a  deeper  more -damag i ng  form  of  fading.  The  correlation 
coefficient  p  =  0  corresponds  to  uncorrelated  (independent)  fading,  while 
p  =  .5  allows  for  adjacent  (equispaced)  pulses  in  figure  1  to  have  some 
degree  of  dependent  fading. 

An  explanation  of  the  initial  result  in  figure  2  follows:  for  K  =  1, 
m  =  l  ,  |_  =  op  (known  noise  level),  the  detection  probability  is  plotted  versus 
the  false  alarm  probability  for  values  of  the  latter  between  IE-10  and  .1. 

Ihe  value  of  the  per-pulse  signal-to-noise  ratio,  E^/N  in  dB,  varies  over 
the  range  6,  8,  10,  ...  ,  42,  giving  detection  probability  values  covering 


*  All  the  figures  are  collected  together  after  the  Summary  section. 
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the  range  .01  to  .990.  The  only  difference  in  the  accompanying  pair,  figures 
3  and  4,  is  that  L  is  reduced  to  32  and  16,  respectively. 

The  results  in  figures  5  through  7  correspond  to  the  worst  cases 
considered  here.  Namely,  there  is  just  one  (fading)  signal  pulse  and  m  is 
.5,  which  means  a  very  deep  fading  medium;  see  [1;  figure  2].  The  values  of 
s igna 1 -to-noi se  ratio  required  for  L  =  16  in  figure  7  are  so  large  as  to  be 
physically  unrealistic,  except  for  the  poorer  quality  region. 

On  the  other  hand,  for  K  =  4  signal  pulses,  Rayleigh  amplitude  fading 
(m  =  1),  and  uncorrelated  fading  (p  =  0),  the  results  in  figures  20  through 
22  are  very  encouraging,  being  physically  reasonable  over  the  whole  range  of 
plotted  values.  But  when  m  is  decreased  to  .5,  and  P  is  increased  to  .5, 
the  results  in  figures  29  through  31,  still  for  K  =  4  pulses,  indicate 
substant ia 1 ly  increased  s igna 1 -to-noi se  ratio  requirements  at  the  higher 
quality  end  of  the  performance  region. 

An  alternative  method  of  presenting  the  graphical  results,  which 

accounts  for  the  losses  incurred  by  not  knowing  the  noise  level,  is  to  plot 

the  required  value  of  E,/N  vs  I.,  for  various  values  of  the  remaining 

1  o 

parameters  and  for  specified  performance  quality  in  terms  of  and  P  . 

Two  such  cases  are  illustrated  in  figures  32  and  33.  They  show  that  the 
cost  of  not  knowing  the  noise  level  is  not  severe  for  the  high  false  alarm 
probab i 1 i t i es ,  but  is  quite  significant  for  the  lower  more-des i rable  false 
alarm  probab i 1 i t i es .  For  example,  in  figure  33  for  K  =  2  signal  pulses,  the 
signal -to -noise  rat  io  must  be  about  1.5  dB  larger  at  L  =  10  noise  pulses  than 
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at  L  =  100,  when  Pp  =  .01.  However,  if  we  want  to  operate  at  Pp  =  IE-10,  the 
increased  signal-to-noise  ratio  requirement  is  about  6  dB  per  pulse.  The 
numbers  are  comparable  for  the  K  =  1  results  in  figure  32. 

The  asymptotes  for  large  L  in  figures  32  and  33  can  be  found  in  some 
cases  from  earlier  results  in  [1],  For  example,  reference  to  [1;  figure  8] 
for  K  =  2,  p  =  .5  gives  E-|/No  -  16.8  dB,  while  Pp  =  IE-6,  Pp  =  .9,  m  =  1. 
Comparison  with  figure  33  here  reveals  that  the  performance  requirement  is 
virtually  at  this  level  by  the  time  that  L  =  100. 
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SUMMARY 


Although  figures  32  and  33  are  very  informative,  allowing  for  a  ready 
assessment  of  the  losses  incurred  by  using  a  finite  small  value  for  L,  the 
number  of  noise-only  pulses,  they  also  illustrate  the  voluminous  compilation 
that  would  be  needed  for  a  thorough  numerical  investigation.  For  example, 
if:  detection  probabilities  PQ  were  of  interest  for  values  .5,  .9,  .99, 

.999;  number  of  signal  pulses  K  for  values  1,  2 . 10;  fading  parameter 

m  for  values  .5,  1,  2;  and  fading  correlation  coefficient  p  for  0,  .5,  1; 
this  would  require  a  total  of  4*10*3*3  =  360  figures.  The  approach  here  is 
instead  to  present  some  representative  receiver  operating  characteri sties,  in 
figures  2  through  31,  from  which  information  similar  to  that  in  figures  32 
and  33  can  be  extracted,  and  to  list  a  general  program  for  the  generation  of 
additional  receiver  operating  characteristics  for  whatever  cases  may  be  of 
interest  to  the  user. 


Some  related  work  on  the  performance  of  a  1  og-norma 1 i zer  subject  to 
Weibul)  or  log-normal  inputs  has  been  published  by  the  author  in  [6]; 
however,  no  fading  was  allowed,  and  the  number  of  signal  pulses  was  limited 
to  K  =  1 .  In  a  different  vein,  the  performance  of  an  or-ing  device 
operating  on  the  output  of  an  incoherent  combiner  of  multiple  pulses  was 
analyzed  in  [7],  These  works  augment  and  complement  the  analysis  conducted 
here . 
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APPENDIX  A 
PROGRAM  LISTINGS 


There  are  two  programs  listed  in  this  appendix,  the  first  for  L  finite, 
the  second  for  L  infinite,  where  L  is  the  number  of  noise-only  pulses  used 
to  establish  a  reference.  The  fundamental  parameters  K,m,l  are  input  in 
lines  20,  30,  40,  while  p  is  input  in  line  1400.  The  particular  values  of 
E^/N  (in  d B )  that  are  of  interest  are  input  in  lines  340  and  350. 

Provision  is  made  for  20  P^  vs  Pp  curves  in  lines  60-90;  this  can  easily 
be  changed  to  accommodate  other  cases. 

The  false  alarm  and  detection  probabilities  are  available  in  lines  1000 
and  1130,  respectively.  The  detection  probability  utilizes  R  and  N  as  input 
variables;  see  table  2.  The  particular  covariance  programmed  in  lines 
1390-1430  is  exponential,  but  this,  too,  can  easily  be  generalized. 

lo  save  space,  the  complete  program  for  L  infinite  is  not  listed. 

Rather,  just  the  essential  false  alarm  and  detection  probability  routines 
are  listed  at  the  end  of  the  appendix;  these  are  obviously  not  functions  of 
L.  The  changes  required  to  accommodate  this  case  of  infinite  L  should  be 
obvious . 
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GENERATE  PD-VS-F'F  NUMBERS  FOR  FINITE  L 

h  =4  '  NUMBER  OF  SIGNAL  PULSES  ADDED,  I 

Ms=.5  1  FADING  PARAMETER,  n,  '.On,  DOF,- 

L  =  1  rf.  i  NUMBER  OF  NOISE  PULSES  ADDED  L 

DIM  U  *  1 0  0  > 

COM  Pf  ■  100  >  ,  Pd  1  1  100',  Pd2 1  100>,  Pd3  <  100>,  P d 4  ■  1  00  >  ,  F  d5  ■  1  00  :• 
COM  PdS <  1  OO  ■ ,  F'd? >  1 OO  ) ,  PdS ',  1  OO  » ,  Pd9 <  1  00  :• ,  Pd  10'  1  00  ■  ,  Pd  11-1  00 
C OM  Pd  1  2  1  00  ■  ,  Pd  1  3  <  1  00  '> ,  Pd  14*'  1  00  >  ,  Pd  1  5  1  1  00  ,•  ,  Pd  1  6  •  1  00  1  ,  Pd  1  7 
COM  Pd  1  3  *  1  00  :■ , Pd  1 9  c  100  >  ,  Pd20  <  1  00  ) 

DOUBLE  K, L, I , J  *  INTEGERS 


FOR  1=1  TO  K 
FOR  J=1  TO  k 
S  =  S  +  F NCo'.'  *  I  ,  J  1 
NEXT  J 
NEXT  I 
k  e  =  y  *  k  S 
N  =  M  =.  *  K  e 
U  =  0. 

U=U+. 01 

Pf  =  FNPt"  *.U,  K,  L  > 

IF  Pf>.l  THEN  200 
U  1  =  M  A  X  C  U  -  .  0  1  ,  .or- 
U  =  U  +  .  O  1 

Pf =FNPt <U,K,L> 

IF  P  f  > 1 E -  1 0  THEN  240 
U2  =  U 

De 1 u=  < U2-U 1 > / 100. 

FOR  I=0  TO  100 
U=Ul+Delu*I 
U(I  '  =U 

Pf ■  I  <  =  F N P f <  U , K , L  ) 

NEXT  I 

FOR  J=1  TO  20 
E  1  r-,odb  =  2>  J  +  2 
E  1  r’i  =10.  '  .  1  *  E  1  ri  o  d  b  > 

R  =  E  1  n  o  *  y  N 
FOR  1=0  TO  100 
U  =  U  •  I  > 

P d  =  F N P d '  U, R, N, \  ,  L  • 

IF  J=1  THEN  Pd  1 '  I  ' =  P  d 
IF  J  =  2  THEN  P  d  2 1  I  1  =  P  d 


THEN  F' d  1  9  1 
THEN  P  d  2  0 1 


1=0  TO  100 
'  =FN  I  rvph  l  ■  Pf  1  I 
I  •  =  FN  I  r  n.'ph  l  P  d  1  ' 
I  >  =FN  I  M"ph  i.  *  F'd2  1 


Pd  19'  I  ■  =FN  I  rv.'ph  l  ■  Pd  19'  I 
Pd20'  I  '  =  F N  I  ri ■  •  p h  l  ■  P d 2 0  1  I 
NE:  T  I 
CALL  A 
END 


'  NORMALIZED  COVARIANCE  COEFFICIENTS 

i  EQUIVALENT  NUMBER  OF  INDEPENDENT  FADES 
!  N  =  m  Ke 


THRESHOLD  VALUES 
P R 0 B A E I L I T V  OF  FALSE  ALAR M 


S I GNAL-T  O-Nu I SE  RATIO  FEP  PULSE,  El  N« 


'  PROBABILITY  OF  DETECTION 


•A 

cr 
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DEF 
IF 

P  =  M  I  N 
T  =  -  L  0  G  <  P  > 

T  =  SuR k T  +  T  ) 

P=  I  . +  T* 1  1.4  3 2 ? 3 
P  =  T-<  2.515517  +  T*' 

IF  X<  .  5  THEN  P=-P 
RETURN  P 
FNEND 

i 

DEF  F  NP  t  IJ  ,  D  0  U  B  L  E  f  ,  L  '  1 

IF  U  =  0 .  THEN  RETURN  1 . 
DOUBLE  Li  1 

U1=U+1  . 
k  1  =  K  -  1 

S=T=EXPkK*LOG<UUl > ) 

FOR  Ls= 1  TO  L-l 
T  =  T*(K1+Ls>^  <  L  i  *  U  1  ) 

3  =  3  +  T 
NEXT  Li 
RETURN  l.-S 
FNEND 

I 

DEF  FNPd  <  U  ,  R  ,  N  ,  DOUBLE  Y  ,  L 
IF  U<=0.  THEN  RETURN  1. 
DOUBLE  Li 
U  1  =  U  +  1  . 

R 1 =R+  1  . 

U  2  =  U  <•’  U  1 
R  u  =  R 1 +  U 

Nk  =N-k  +  1 
Y  =  R  1  R  u 
X=U2*R' R 1 
X  1  =  X  -  1  . 

:  =  T - E X P 1  K ♦LOO1  U 
Ho  =  U  . 

H=  1  . 

FOR  L i = 1  TO  L-l 
T  =  T*Y 
.1  =  k.  2  +  L  i 

H=  <  •:  J  +  L  i  +  1  N(  -Li 
Ho  =  H 
I  H  =  Pl 
i  'j  =  b  +  T  ♦  H 

NEXT  Li 
RETURN  1.-3 
FNEND 


i  t  T  t ,  ij  ki  l  3  ij  >: 

0  1  0  3  2  i  •  ■  P 

F H L  ■■  E  hlhF  M  P R  0  B Pi B  I  L  I  T  i 

INTEGER 


1  DETECTION  PROBABILITY 
i  INTEGER 


1  +  N  *  L  0  G  k  U  1  R  u  >  > 


)*h+ j*x.i*ho *  Li 


F  N  I  n  k1  ph  l  1  X  1  1  H  M  3 

5  THEN  RETURN  0. 

X ,  1  .  - .  1 


+  T ♦ ■  .  1392 6 
O  2  3  5  3  +  T 


DEF  FNC c"  DOUBLE  I  ,  J  ■ 

R h =  .  5  1  NORMALIZED  COVARIANCE  COEFFIC  IENT 

C  o'<  =  F’ho  '  HB3  ■  I  -  J  1  '  EXPONENT  I  AL  BEHAVIOR 

R'  E  T  U  R  N  C  o 1 1 
FNEND 
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•a 


SUB  fl  1  PLOT  PD  VS  PF  ON  NORMAL  PROBABILITY  PAPER 
C  0 M  P f 1  *  ; ,  P d  1  * *  ; ,  P d 2  *.  *  ; ,  P d 3  *  *  * ,  P  d  4  *  *  > ,  P d  5  <  *  • 

COM  P d 6  *  *  ;• ,  Pd? *.  *>,  PdS *  * ; ,  Pd9 *  *  > ,  Pd  1 0 * .  *  > ,  Pd  1  1 1  *  • 

COM  Pd  1  2  <  *  ' ,  F'd  1  3  v  *  > ,  Pd  1  4  <  *  ;  ,  Pd  1  5  *  *  •  ,  Pd  1  6  ■-  -  '  ,  P  d  1  ?  * .  *  - 
COM  PdlS* »  1 , Pdl9( *  > , Pd2Q(* ) 

DIM  Pit  30] , B* C  303 
D  I  M  X  label  1 <  1  :  3  0  ;  ,VUbel*vl:  3  0  ) 

D I M  Xc oord 1  1  :  30  > ,  Yc oor d* 1 : 30  » 

D  I  M  X  g  r  l  d  *  1  :  3  0  ;  ,  Y  g  r  l  d  *  1:30  1 

DOUBLE  N ,  L  x  ,  L  y  ,  N  x  ,  N  y  ,  I  !  TWTES-ERS 

i 

A* = " Probab ll i t y  of  False  Alarm" 

B *  =  " P r o b ab i  1  i r y  o f  Detectio n " 

i 

L  x  =  1  2 

P  E  D I M  X  label**  1 : L x  > , X c o o  r  d  t 1 : L  , > 

DATA  E-10,E-9,E-8,E-?,E-6,E-5,E-4,E-3, .01, .02, .05, . 1 
READ  X  label***; 

DATA  IE-10,  IE-9,  IE-8,  IE-7,  IE-6,  IE-5,  IE-4,  .001,  .01,  ,02,  .05,  .  1 
READ  Xc  oor d ( * ) 

i 

L  y  =  1  a 

R  E  D  I  M  Y  1  ab  e  1  *  <  1  :  L  y  >  ,  Y  c  o  o  r  d  *  1  :  L  y  > 

DATA  . 0 1 , . 02 , . 05 , . 1 , . 2 , . 3 , . 4 , . 5 , . 6 , . ? , . S , . 9 
DATA  . 35 , . 98, . 99, . 995 , . 998 , . 999 
READ  Y label***; 

DATA  . 0 1 , . 02 , . 05 , . 1 , . 2 , . 3 , . 4 , . 5 , . 6 , . ? , . 3 , . 9 
DATA  .95, .98, .99, .995, .998, .999 
READ  Y coord**; 

i 

N  x  =  1  4 

REDIM  X grid* 1 : N  < ; 

DATA  IE-10,  IE- 9,  IE-6,  IE-7,  IE-8,  IE-5,  IE-4 
DATA  . O0 1 , . 002 ,  . 005 ,  . 0 1 ,  . 02 ,  . 05 ,  .  1 
READ  1 !  grid* * ; 

i 

Ny=  1  8 

R  E  D  I  M  Y  g  r  i  d  ■  1  :  N 1  >  > 

DATA  .  0  1  ,  .  02  ,  .  05  ,  .  1  ,  .  2  ,  .  3  ,  .  4  ,  .  5  ,  .  6  ,  .  ?  ,  .  8  ,  .  9 
DATA  . 95, . 98, . 99, . 995, . 998, . 999 
READ  Ygr i d 1  *  > 

i 

FOR  1=1  TO  Lx 

1  oor  d  *.  I  ;  =FN  I  rv'ph  i  <  Xc  oor  d  *  I  1  1 
NEXT  I 

FOR  1=1  TO  Ly 

Y c  o o r  d  *  I  ;  =  F  N  I  n  v p  1;  t  *  Y c  o or  d  *  I  >  1 
NEXT  I 

FOR  1=1  TO  Nx 

,'gr  id1  I  ■  =FN  I  rtvpb  i  1  Xgr  id1  I  1  > 

NEXT  I 

FOR  1=1  TO  Ny 

Ygr id*  I  ;  =FN I nyph l  *  Ygr id'  I  >  ■ 

NEXT  I 
X  1  =Xgr i d  *  1  1 
X2  =  Xgr id*  U  -  > 

Y  1  = Ygr i d'  1  • 

*'2  =  'i'gr  id1  N"  * 

e  =  ■  Y  2  -  v  l  ■  :  2  - :  l  • 
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GINIT  200. 

2  E  0 . 

'  VERTICAL  PAPER 

PLOTTER  13 

505.  "  HPGL " 

PRINTER  13 

505 

PRINT  " V  S  4 " 

LIMIT  PLOTTER  5  0  5 , 0  .  ,  2  0  O  .  ,  0  .  ,  2  6  0  .  1  1  G  B  U  =  2  i  n  in 

VIEWPORT  22  .  ,  85.  ,  13.  ,  122. 

VIEWPORT  22. ,35. ,59. , 122.  1  TOP  OF  FhPER 

VIEWPORT  22.  .  35.  ,  19.  ,  62.  1  BOTTOM  OF  F'FiPEP 

WINDOW  X 1 , X 2 , Y 1 , V 2 

FOR  1=1  TO  Nx 

MOVE  Xgridi  I  >,  Y  1 

BRtiW  Xgr i d 1  I  '  ,  Y2 

NEXT  I 

FOR  1=1  TO  Ny 
M  0  V  E  XI, Y  gr i d (  I  > 

DRhW  X 2 , Y g r t  d < I > 

NEXT  I 

CSI2E  2. 3,  .  5 
LORG  5 

V  =  Y  1  -  (.  Y2-Y  1  >  *  .  02 

FOR  1=1  TO  Lx 
MOVE  Xc  cor'd1-  I  1  ,  Y 

LHBEL  X 1 db e 1  $ i  I  ' 

NEXT  I 
CSIZE  3.  ,  .  5 

MOVE  .  5*  <.  X  1  +X2  }  ,  Y  1  -  .  OE*  ‘  Y2-Y  1  > 

LHBEL  H* 

MOVE  ,5*i  X 1  +X2  > ,  Y 1  -  .  1  *  <  Y2-Y  1  '> 

LHBEL  "Figure  31.  ROC  for  K  =  4,  in  =.5,  =.5,  L  =  1  E " 

CSI2E  2. 3, . 5 
LORG  S 

::  =  :< l - >' X2-X l  '*.01 
FOR  1  =  1  TO  L y 
MOVE  X , Yc oor  d  >  I  > 

LhEEL  Y 1  abe 1  *  I  • 

NEXT  I 
LBIR  PI  2. 

CSIZE  3. , . 5 
LORG  5 

M  0  V  E  X1-.  15*'  :  -  - .  1 1  .'  ,  .  5  *  '-  Y  1  +  Y  2  1 

LHBEL  Bt 

PENUP 

PLOT  P  t  i *  * , F  d  1  '  *  ' 

PENUP 

PLOT  P  t •  *  < , P  d  2 i *  > 

PENUP 


PLOT  P  t  >  *  j  ,  P  d  1  ’f"  *  ■ 

PENUP 

PLOT  P  t'  i  *  :>  ,  P  d  2  0  1  *  > 

PENUP 

PhU'GE 

PRINTER  13  CRT 

plotter  505  is  terminated 

3 UBEND 
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